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1 Introduction 



In the study of self-gravitating systems there are three factors whose relevance 
has been recurrently and separately stressed in the literature. These are, the 
Weyl tensor, the local anisotropy of pressure and the inhomogeneity of energy 
density distribution (density contrast). 

The Weyl tensor |L| or some functions of it [0, have been proposed to 
provide a gravitational arrow of time. The rationale behind this idea being 
that tidal forces tend to make the gravitating fluid more inhomogeneous as 
the evolution proceeds, thereby indicating the sense of time. However, some 
works have thrown doubts on this proposal || . Also, as it will be seen below, 
it is worth noticing that the relation between the Weyl tensor and the density 
contrast is affected by the presence of local anisotropy of pressure. 

The role of density inhomogeneities in the collapse of dust 0] and in 
particular in the formation of naked singularities ||, has been extensively 
discussed in the literature. 

Finally, the assumption of local anisotropy of pressure, has been proved 
to be very useful in the study of relativistic compact objects (see |5| and 
references therein). 

A hint pointing to the relevance of the above mentioned three factors in 
the fate of spherical collapse is also provided by the expression of the active 
gravitational mass in terms of those factors [0], f|. 

These three factors are usually considered separately, their relationship 
being omitted from discussion,even though they are related by a simple ex- 
pression, which we shall present below M, ||. 

The purpose of this work is twofold. On one hand we shall integrate the 
vanishing Weyl tensor condition, which will allow us to construct conformally 
flat models (with anisotropic pressure). The obtained solutions represent 
static or slowly evolving (in the quasi-static approximation) spheres, which 
could serve for the modelling of compact self-gravitating objects. On the 
other hand, we want to study, comparatively, the effects of the above men- 
tioned parameters on the (slow) evolution of relativistic spheres. For doing 
this we shall contrast the evolution of the conformally flat models with the 
evolution of models with the same energy density or radial pressure distribu- 
tion, but non-vanishing Weyl tensor. With this purpose it will be useful to 
calculate the active gravitational mass and the fluid velocity for each model. 

The paper is organized as follows: In the next Section all relevant equa- 
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tions and conventions are given. The condition for the vanishing of the Weyl 
tensor is integrated in Section 3 and the models are described in Section 4. 
Finally a discussion of results is presented in the last Section. 



2 Relevant equations and conventions 
2.1 The field equations 

We consider a spherically symmetric distribution of collapsing fluid, which 
we assume to be locally anisotropic and bounded by a spherical surface E. 
The line element is given in Schwarzschild-like coordinates by 

ds 1 = e u dt 2 - e x dr 2 - r 2 (d9 2 + sin 2 9d(t) 2 ) (1) 

where v and A are functions of t and r. The coordinates are: x° = t; x 1 = 
r; x 2 = 8; x 3 = 0. 

The metric ([!]) has to satisfy Einstein field equations 



G$ = -8tKC (2) 



which in our case read 



^ = --,+z- x \ir,-'-\ <■>> 
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SnT 2 = -8ttT| = - — (2A + A(A 



^ ( 2v" + v' 2 - AV + 2^^) (5) 



-87tT i = -- (6) 
r 

where dots and primes stand for partial differentiation with respect to t and 
r respectively. 
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In order to give physical significance to the components we apply the 
Bondi approach 0, i.e we introduce local Minkowski coordinates (r,x,y,z), 
defined by 



dr = e u/2 dt dx = e x/2 dr 



dy = rdO dz = rsin 



Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



rpO 

1 o 



Tl = T\ 



Tl 



T 



01 



-{v+X)/2 T 



01 



Next we suppose that, when viewed by an observer moving relative to these 
coordinates with velocity u in the radial direction, the physical content of 
space consists of an anisotropic fluid of energy density p, radial pressure P r 
and tangential pressure P±. Thus, when viewed by this moving observer, the 
covariant energy-momentum tensor in Minkowski coordinates is 



(P 

P r 

Pi 

V o o o 



\ 




Px J 



Then a Lorentz transformation readily shows that 

p + P r u 2 



J n — -Ifl 



Tl = Tl 



Tl = Tl = Tl 



l-uo 2 

Pr + pU 7 
1 - 

= f, 3 



0J* 



Pi 



T 01 = e^ 2 f 01 : 
Note that the velocity in the (t, r, 9, 
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dr 
dt 



{p + P r )ue^ +X ^ 2 
1-uj 2 

system, dr/dt, is related to uj by 

,(\-v)/2 



(7) 

(8) 
(9) 
(10) 

(11) 



Outside of the fluid, the spacetime is Schwarzschild, 



ds 2 



2M\ 



dt 2 - 1 



2M 



dr 2 - r 2 (d6 2 + sin 2 6d<f) 2 ) (12) 



In order to match the two metrics smoothly on the boundary surface r = 
rs(t), we require continuity of the first and second fundamental forms across 
that surface. As result of this matching we obtain the well known result 

[Pr\x = (13) 

Next, the radial component of the conservation law 

= (14) 

gives 

167T 



{-&xTi 




Tl-Tl +W [Tl-Tl ) + — X + -ir-Z) (15) 



which in the static case becomes 

^ / \ 2 (P ± - P r ) . . 

K = ~^(P + Pr) + ^ \ ~ (16) 

representing the generalization of the Tolman-Oppenheimer-Volkof equation 
for anisotropic fluids ||. 

In this work we shall consider exclusively static or slowly evolving (quasi- 
static) systems. By this we mean that our sphere either does not change or 
changes slowly on a time scale that is very long compared to the typical time 
in which the sphere reacts to a slight perturbation of hydrostatic equilibrium 
(this typical time scale is called hydrostatic time scale). Thus our system is 
always very close to or in hydrostatic equilibrium and its evolution may be 
regarded as a sequence of static models linked by (||). This assumption is 
very sensible because the hydrostatic time scale is very small for almost any 
phase of the life of the star. It is of the order of 27 minutes for the Sun, 4.5 
seconds for a white dwarf and 10~ 4 seconds for a neutron star of one solar 
mass and 10 Km radius. It is well known that any of the stellar configurations 
mentioned above, change on a time scale that is very long compared to their 
respective hydrostatic time scales. Let us now translate this assumption in 
conditions to uj and metric functions. 

First of all, slow contraction means that the radial velocity u as measured 
by the Minkowski observer is always much smaller than the velocity of light 
(w < 1). Therefore we shall neglect terms of the order O(uo 2 ). 
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Then ( |1"5|) yields 



v A 2 z>A 

A + = 

2 2 



+ (p + P r ) - - 2 



v' Pi -P 



2 r 



(17) 



Since by assumption, our system is always (not only at a given time t) 
in equilibrium (or very close to), ( |16|) and fll7|) imply then, for an arbitrary 
slowly evolving configuration 

A « z>A « A 2 « (18) 

and of course, time derivatives of any order of the left hand side of the 
hydrostatic equilibrium equation must also vanish, for otherwise the system 
will deviate from equilibrium. This condition implies, in particular, that we 
must demand 

v m 

Finally, from the time derivative of @, and using (0), it follows that 

Co « 0(A, \lu,vlu) (19) 

which implies that we shall also neglect terms linear in the acceleration. On 
purely physical considerations, it is obvious that the vanishing of Co is required 
to keep the system always in equilibrium. 

Thus, from now on, we shall always assume 

0(u 2 ) = \ 2 = u 2 = \u = \ = i) = (20) 

implying that the system remains in (or very close to) equilibrium. 



2.2 The Weyl tensor 

We can now calculate the components of the Weyl tensor. Neglecting terms 
of order Az>, A 2 , z> 2 , A and i), we find that all non-vanishing components can 
be expressed through C| 32 . Thus, 

r , r 3 e" A fe x v'\' 1 v' 2 v" \' - v'\ . . 
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Next, defining the mass function as usual 

m(r,t) =4tt f T°r 2 dr (22) 
Jo 



the following relations may be established (0],Q 



w 



-l*lr*(n)'dr + l<Kr*(Tl-Tl) (23) 

m(r,t) = ^vrr 3 T ° - ^ jTr 3 (T °)'dr (24) 

Both, (p3|) and ([24]) are valid in the general (dynamic) case. However only 
in the static or the quasi-static case Tq and T\ denote the proper energy 
density and the radial pressure respectively 

2.3 The Tolman-Whittaker mass 

The Tolman-Whittaker mass JTO] within a sphere of radius r inside S, is 
defined as || 

miw (r, t ) = 4tt r 2 e^ +A )/ 2 (t o ° - - 2T 2 2 ) dr (25) 

Two alternative expressions, easily obtained from the field equations (see 
[§] for details) are 

m TW = (>+ A )/ 2 + 47rP r r 3 ) (26) 

and 

TfiTw = ie<"- A >/Vr 2 (27) 

The interpretation of mrvy as the active gravitational mass follows at 
once from (p7 ) and (p!6|). Indeed, the first term on the right side of (p~6|) (the 



"gravitational force" term) is the product of the passive gravitational mass 
density (p + P r ) and a term proportional to rriTw/ r2 - A similar conclusion 
can be obtained if we recall that the gravitational acceleration of a test 
particle, instantaneously at rest in a static gravitational field, as measured 
with standard rods and coordinate clocks is given by [O 



eC- A >/V m TW 

a = = — 28 

2 r 2 K ' 
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2.4 The velocity of a fluid element 

For the comparative study of the (slow) evolution of different solutions, it 
will be useful to plot the velocity (u) profiles for different pieces of material. 
A simple expression for uj, may be obtained as follows: from (^) and (P%fl, it 
results 

(29) 



—A 



2m 



and from (|29|) 



then using (|6]) 



Are" 



m 



me 



(X-v)/2 



UJ 



Anr 2 (p + P r 



(30) 



(31) 



3 The vanishing Weyl condition 

We shall now proceed to integrate the condition 

W = 



which, using (21), may be written as 

+ e 



e- x - 1 



2r 



Introducing new variables 



2r 



y = e 



equation(^) is cast into 



^ , + 2y [u" - u'/r + u/r 2 ] 



~2 ~ ~u 



2u 



[u' — u/r] 



r 2 [u' — u/r] 



whose formal solution is 



y = e 



— J k(r)dr 



zS k ^ dr f{r)dr + d 



(32) 



(33) 



(34) 



(35) 



(36) 
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where G\ is a constant of integration, and 



, / % d 
k(r) = 2— 
dr 



f(r) 



in I u 

V r 

2u 



r 2 (u' — u/r) 

changing back to the original variables, eq.(|36|) becomes 

u> 1 p A / 2 



c 2 r 2 e v 



with c 2 = —C\. 

Next, (|39|) may be formally integrated, to obtain 



c 2 r 2 cosh^ 



3 A/2 



-dr + C 



(37) 
(38) 

(39) 



(40) 



where C is a constant of integration (a function of t, in the slowly evolving 



case). The reader may check that ([40]) satisfies ( p9|) and (|32D (or (33 )). 
In the next section we shall present some models satisfying (R^) (or (pEUD 



4 The models 

A simple counting of equations ((H)-©) and unknowns (u, A, p, P r , P±, uj) 
indicates that we have to provide two additional relations (in the form of 
equations of state and/or restrictions on metric variables), in order to inte- 
grate the system (@ _ (H))- 

If one assumes that the fluid is locally isotropic (P r = P±) then, demand- 
ing W = 0, we are driven to a unique solution (the Schwarzschild interior 
solution), a fact also obvious from (|2"3"D. However, if P r ^ P±, then the 
condition W = 0, does not single out a unique model. 

In what follows we shall construct two models with W = 0. One of them 
characterized by P r = 0, and for the other we prescribe a given energy density 
distribution which is similar to the one proposed by Gokhroo and Mehra [1^ 



Additionally, we present two other models with W ^ 0. One is characterized 
by P r = and p = p(t) (|P|,[^3|), and the other has the same energy density 



distribution as one of the conformally flat solutions. 
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4.1 Model I 

Our first model is characterized by 



W = 



and 



P r = 

Then, from (H) and (0), it follows 



(4-9</) 



with 



c 2 r 2 _ 2p ) 
__ m(r, t) 



and where (|29|) has been used. 

Next, taking r- derivative of fl43"D , and using 



2m 



20 



r 2 (l — 2m/r) r(l — 2g) 
easily derived from ([|), we obtain 

540 3 - 420 2 + 8c/ 







which after integration yields 



r (I80 2 - 180 + 4) 



Dr 



.1/2 



(4-90)V6 



(41) 
(42) 

(43) 
(44) 



(45) 



(46) 



(47) 



where D is a constant (a function of t in the slowly evolving case) of integra- 
tion. Solving (|47D for g, one obtains 



= a 



1/3 



with 



2 + (27a + 4) 1/2 
(Dr) 6 = a=( 



1/3 



+ 



2 - (27a + 4) 



1/2' 



1/3' 



9l 



r^J 4 -90s 



(48) 
(49) 
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where subscript X indicates that the quantity is evaluated at the boundary 
surface r = rs- 

The remaining variables are now easily obtained from the field equations 
and (H), (H). Thus, 

3£ (1-20) 



P 
P, 



2vrr 2 (2 - 3#) 
3^ 2 1 



4vrr 2 (2 - 3#) 
e- X = 1 - 2 5 

For the Tolman-Whittaker mass we obtain, using either (|26|) or 

1/2 



m T1 y = # 



9e 
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Is) 



(50) 

(51) 
(52) 

(53) 



or, using the dimensionless varibles 

r 



x = 



m s M 
n = — = — = 



Mi 3 (Zi + Z 2 ) 1/2 (4-9n) 
4(4-9n) 1/2 -9m; 2 (Z 1 + Z 2 



,1/2 



with 



2 (4 - 9n) 1/2 ± (27xV + 16 - 36n 



1/2 



1/2' 



1/3 



(54) 
(55) 

(56) 



It is worth noticing that from the requirement p > 0, it follows, using 



9< 3 



a stronger restriction appears from the condition 

P>P± 

which requires 

2 

q < - 
y 5 



(57) 

(58) 
(59) 
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Finally, the velocity 00 of any fluid element is given by 
u^n 1 ' 2 x 2 {2 - 3g)(4 - 9g) 1/2 (l - 2n) 

^ 1 /2(4_ 9n )2( 1 _ 2 ^)2 X 

x {2(Z 1 + Z 2 ) + ^ ± -m\ 60 

I (27x 6 n 3 + 16-36n) 1/2 J 

where (EH), fBl, fEl)-(l50|), f|5l, and (M) have been used. 



4.2 Model II 

The second model we shall consider, is well known. It is characterized by 

P r = (61) 
P = PH) (62) 

and is not conformally flat. 

A static version of this model was studied by Florides p~3|. Models with 



vanishing radial pressure have been discussed in the past |fL4| , and more 
recently in relation with the formation of naked singularities [To* . 
The corresponding variables are (see for details) 

p = p(t) ; P r = (63) 

r,-^^ (64) 



3 (l - f r 2 p) 



1 ~ f rip 
1 - f r 2 p 



3/2 



1/2 



(65) 



e- x = 1 - ^r 2 p (66) 
ir (67) 



9(1- ^r 2 p 



47T 



8^2 „\ 3/4 



—-f^(MS) (68) 



/v / \ 1 - ^r 2 p 



- = I^-*St] (69) 



12 



4.3 Model III 

This model is conformally flat and is further characterized by 

/ „2\ 2 



1 



b 2 



(70) 



where b is a constant (a function of t in the slowly evolving case). 
Then, from (EDI) one obtains 



c 2 [B 2 r 2 + b 2 -r 
4B 2 (b 2 - r 2 ) 



21 2 



with 



C = InB 



(71) 



(72) 



^From the field equations, ([71]) and ([H]), we can now obtain the expres- 
sions for p, P r and P±. 



P 



Anb 2 



5r 2 



SnP r = -— + 



1 (b 2 - r 



2\2 



1 AB 2 b 2 -2[B 2 r 2 + b 2 -r 2 



b 4 

8tt {P r - P 



+ 



^■•2 ^^2 ^ 2^2 j 7™ J 

_2r2 

6, c and £? are related to the total mass and the radius of the sphere through 
the boundary conditions 



(73) 
(74) 
(75) 



1 - 



2M 
2M 



P rE =0 



the corresponding expressions are 



2M 



1/21 V2 



(76) 

(77) 
(78) 

(79) 
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AM 9M 2 



1/2 



5 



;i - 2M/r s ) 1/4 [(1 - 2M/r s ) 1/2 + (3M/r E - 1) 



1 - (1 - 2M/r s ) 1/2 j 1/2 [4M/r E -9M 2 /r|] 1/2 



(80) 
(81) 



For the Tolman-Whittaker mass the obtained expression is 

„3 



cr 



rriTw 



2Bb 2 {b 2 - r 2 ) 1 / 2 
or, using ©-(HD and (H 

Mx 3 (l - 2nfl i 



2B 2 b 2 - B 2 r 2 -b 2 + r 2 



(82) 



m TW 



3-2x 2 



[1 -x 2 + x 2 (l -2n) l / 2 ] 1/2 
(1-x 2 



x 



n 



1 - (1 - 2n) 



1/2 



1 + 



An — 9n 2 



2(l-2n)V2 [(i _ 2n) 1 /2 + 3n _ i] 
and the expression for uo in this model, results in 



AbBx 



1 - (1 - 2n) 



1/2 



(An - 9n 2 ) 1/2 (1 - x 2 [1 - (1 - 2 n y/ 2 }) 1/2 

1 

X [25 2 + B 2 x 2 [1 - (1 - 2n)V2] + 1 - [1 - (1 - 2n) 1 /2] 
where 6 is easily obtained from ( [79] ) 

. u; s n(l - 2n)V2 (5 — ^ [l — (1 — 2n) 1 / 2 ] ) 



(84) 



2 [1 - (1 - 2ny/ 2 f 2 



(85) 



4.4 Model IV 



This model has the same energy density distribution as the previous one 
(same A), but is not conformally flat (W ^ 0). Instead, the model is further 
characterized by 

P r = (86) 
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Then, from flT0|), ( p6[ ) and field equations, we obtain 

3 / 5r 2 \ 



47Tb 2 



6b 2 I 



, 2(b' I -r' 1 ) 



with 



(2b 2 r 2 - r 4 ) (6fe 2 - 5r 2 ) 
32tt6 4 (6 2 - r 2 ) 2 

2M\ 5 / 4 



2(l-2A//r E ) 1 / 2 



(87) 

(88) 
(89) 

(90) 



and b is given by (|79|). 

For the active gravitational mass, we obtain, after some lengthly calcula- 
tions 



m TW = M(l - 2n) 5/ V (l - (1 - 2n) 1/2 



1 _ (l _ (1 _ 2n) 1 ' 2 



x exp - 
4 



(l-x 2 [l-(l-2n)V2]) 5 / 4 
(a; 2 - 1) (l - (1 -2n) 1 / 2 ) 



(1 - 2n)V2 [l _ X 2 (i _ (i _ 2n)V2)] 



(91) 



where n and x are defined by eq. (|54|) . Observe that condition p > is 
satisfied for all n, however if we demand p > P±, then 

n < 0.4 (92) 

Finally, the expression for the velocity takes the form 



CO 



2uj^x l-x 2 + x 2 {l-2nf' 2 5n - 2 + 2(1 - 2nfl 2 



1/4 r 



(1 - 2n) 1 /8 [ 6 
(1 



cxp ■ 



5a; 2 + 5x 2 (l - 2n) 1 /2] [1 - (1 - 2n) 1 /2] 
x 2 ) [l-(l-2n) 1 /2 



(1 - 2n)V2 [l-x 2 + x 2 {l - 2n)V2] 



(93) 
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5 Discussion 



We have integrated the vanishing Weyl condition. The resulting expres- 
sion (|39| ) (and (|40|) ) allows to find conformally flat models in a very simple 
way, once an additional condition on physical or metric variables is imposed. 
Specifically we have found two conformally flat models (I, III). In order to 
bring out the role of Weyl tensor in its slow evolution, we have also pre- 
sented two other models (II, IV) with non-vanishing Weyl tensor. Model II, 
as model I has vanishing radial pressure, whereas model IV has the same 
energy density distribution as model III. This will allow to see the effect 
of the three abovementioned factors (local anisotropy, Weyl tensor, density 
contrast) on the Tolman-Whittaker mass distribution within the sphere, and 
on the velocity profile of different pieces of matter. With this purpose, all 
models are considered with the same total mass M and surface velocity ws. 

Figure 1 exhibits the evolution of ttitw/M as function of x, in the process 
of slow contraction (increasing n), for the model I. For all other models the 
behaviour is qualitatively the same, i.e. as the contraction proceeds, the 
active gravitational mass within the sphere decreases. However, the absolute 
value of rriTw is different for different models (for same n and x) as can be 
seen from figures 2-4, which display the ratio vtitw / /m TwO-ty f°r the three 
models (I, III, IV). 

As it can be seen, for any r < (for the same total mass M), we have 

mTiy(IH) > m TW (l) > m TW (IV) > m TW (II) 

the differences being larger for more compact (larger n) configurations. Par- 
enthetically, the two conformally flat models present the largest TW masses. 

For models II, III and IV the collapse proceeds in a quasi-homologous 
(quasi-linear) regime as indicated in figure 5 for model III (for models II and 
IV the figures are similar), deviating from that regime as n increases. 

However for model I, the contraction is not homologous even for small n 
as indicated in figure 6. 

Figures 7-9 display the ratio u>/u(II) for the three models (I, III, IV). 

Except for extremely high fields in model I, we see that 

cj(IV) > w(III) > w(H) > w(I) 
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which indicate that for the same energy density distribution (III, IV) or radial 
pressure distribution (I, II), the slow contraction of interior shells proceeds 
slower in the conformally flat case. 
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6 Figure captions 

Figure 1: ititw/M as function of x for the model I, and ten values of n, from 
.04 to .4. 

Figure 2: m TW (I) /m TW (I I) as function of x , curves a — j correspond to 
n=A; .36; .32; .28; .24; .20; .16; .12, .08; .04. 

Figure 3:mTw(IH)/ m Tw(H) as function of x ,for the same values of n as in 
figure 2. 

Figure 4: mTw(IV) / m Tw(I I) as function of x ,for the same values of n as 
in figure 2. 

Figure 5: u/uj^ as function of x for the model III, and ten values of n, from.04 
to .4 

Figure 6: Same as figure 5 for the model I. 

Figure 7: The ratio cu(I)/cc;(II) as function of x, for ten values of n from .04 
to .4. 

Figure 8: The ratio uj(1V)/uj(1T) as function of x, curves a — j correspond to 
n=.4; .36; .32; .28; .24; .20; .16; .12, .08; .04. 

Figure 9: The ratio u^(III) / u^(II) as function of x, for ten values of n from .04 
to .4. 
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